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ABSTRACT 

This paper considers problems which accompany thermionic 
emission of electrons from a hot body surrounded by a plasma. 
In the absence of other mechanisms, an electric potential is 
established at the surface of the body through the balance of. 
thermionic emission and accretion of electrons from the 
external plasma. Analytical solutions are obtained for the 
electric potential field and the electron density distribution 
around the body. A possible application of this analysis to 
objects in space is indicated. 
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I . INTRODUCTION 

An object in space may become hot while approaching a 

hot stellar body like the sun, or while entering a dense 

atmosphere like that of the earth. Long before such metallic 

objects melt, evaporate, or ablate, they may acquire temperatures 

which are sufficient to cause a copious emission of electrons 

from the surface. Therefore, the temperatures . lower than, and 

in the neighborhood of, the melting point are of interest to 

us in this paper. As a matter of convenience and without serious 
y 

loss of generality, we will regard iron as a reference' substance 
composing the objects in space, and hence consider temperatures 
lower than 1600°K. The analytical formulae are, however, 
applicable to any other specific case of a surface capable of 
thermionic emission. 

Thermionic emission is very sensitive to temperatures; the 

- : ~ ' : Vo " lg' 

emitted electron flux is of the order of 10 and 10 

electrons/cm^-sec at surface temperatures of 1000°K and 1600°K 
respectively / from a material of work function W_ = 3 ev. The 

v u 

#' 

emission of electrons from the object's surface leaves a positive 

surface charge. A great majority of the emitted electrons 

\ 

describe ballistic orbits and return to the surface, while a 
certain number of those in the high energy tail of the energy 
spectrum are able to escape from the potential field of the 
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object. The positive charge at the object's surface is 
established by these escaping electrons, and the rate of etecape 
of tho emitted electrons decreases with an increase in the 
surface potential. Furthermore, if the object is surrounded by 
a plasma, the plasma accretion alone has a tendency to impart 
a negative charge to the object's surface. Therefore, a steady 
potential can be established at the object's surface when the 
net negative charge leaving the object due to thermionic emission 
is completely replenished by the net negative charge brought to 
the surface by the accretion from the surrounding plasma. The 
magnitude of the equilibrium surface potential is then determined 
from the balance of the plasma accretion current and the escape 
component of the thermionic emission current. 

There are other mechanisms (Chopra 1961) in which an object 
may acquire an electric charge. An effect of considerable interest 
is connected with the photoelectric emission and accretion of 
electrons. The photoelectric effect is important for objects on 
the day side of the earth and for surfaces exposed to the sun. 

In certain cases, it is comparable to, and at times may even become 

V 

more significant than, the thermionic emission. We will, however, 
limit the analysis of the present paper to only thermionic emission 
and leave these other considerations for a subsequent paper. 

The incoming plasma electron -flux and the thermionic electrons 


* 
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constitute a plasma cloud with most of the contribution to the 
electron density in the cloud coming from the ballistic component 
of the emitted electrons. This plasma cloud screens the electric 
potential on the body. An analytical expression for the density 
distribution in terms of potential cp 0 and work function W Q is 
■ obtained by solving the equations of Poisson and the conservation 

of energy and momentum. This analytical expression is substituted 
back in the Poisson equation which is then solved numerically to 

! yield potential distribution as a function of distance from the 

i 

surface. 

| r* 9 

l It may be mentioned here that the problem considered in thxs 

| paper bears a certain analogy to the problem of the exosphere* 

i 

j In the exosphere problem, the particles are projected outwards 

| corresponding to the temperature of the base layer. One of us 

(C.S. Shen) has successfully applied (Shen 1963) the present 
• analysis (after some modifications) to the structure of the 

j planetary exosphere, and has obtained an analytical expression 

for the density distribution. 

V 

II. FORMULATION OF THE PROBLEM - BASIC EQUATIONS 

it 

Let us consider a spherical object with an equilibrium 
surface potential cp and surface temperature T, surrounded by 


y 
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screening charges due to thermionic emission and a rarefied 
external plasma with electron density n^ and ion density n Q /z 
(where 2 e is the ionic charge at a temperature T p X When the 
( thermionic emission is stronger than the plasma accretion and 
the object is moving slower than the mean thermal speed of the 
plasma 'electrons (~10 7 cm/sec) the~>otent£al~co{f ) and the 
screening electron density p{r) are, to a first degree of 
approximation, spherically symmetrical, and are given by 

7 3 cp(r) = - f p(r) U) 

where e = 4.8 x 10 _1 ° e.s.u. is the electron charge, = 1 
is the permitivity of the medium, and r is the radial distance 
measured from the center of the spherical body. 

The screening electron density p(r) consists of three parts 

p(r) = p b (r) + P esc (r) + p p (r) . (2) 

Here p, (r) is the ballistic component which is due to the electrons 
emitted from the surface with velocities less than the escape 
velocity; these particles describe ballistic orbits in the electric 
potential field of the body and return to the surface. The escape 
component p (r) is due to the electrons emitted with velocities 

0SC 
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exceeding the escape velocity; these particles do not return 

to the charged body. The third component p(r) is due to the 

accretion from the surrounding plasma. Among these (r) 

contributes about SO percent to the local electron density (as 

can be seen from later calculations) . Also, in the steady state 

condition, the escape component of the thermionic electron flux 

is equal to the incoming plasma accretion flux. Therefore, to 

simplify one of our later calculations, we can set o (r) = n (r) 

p esc K p 

Assuming that the electrons inside the metallic body have 


velocities given by the Fermi distribution law, the number of 
electrons having velocities in the range (v, v'+ dv) and hitting 
a unit area of the surface (inside) is given by 


4 rim 3 v v dv dv 
J (y) = — r 


(E-E f ) AT +1 


(3) 


where 


E = l/2m(v 2 + v 2 ) 
r t 


, -28 -27 
and m = 9 x 10 g is the electron mass, h = 6.27 x 10 erg/sec 

is the Planck's constant, E f is the Fermi energy, and v r and v 

are the components of the velocity v in directions parallel and 

transverse to the radius vector r. 


If we denote the velocity of the electron at the position 
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r, (r>R) , by u (r) , then the principles of conservation of energy 

' (4) 

(5) 


( 6 ) 

which provides a stringent condition for an electron emitted = 
from the surface to reach the radial distance r. Only those 
electrons with initial velocity v. satisfying the inequality 

v r 2 + (1 - | 2 >v t 2 - I f e (cr> 0 - oj + E f + W Q } > 0 (7) 

can reach position r. These electrons may be divided into 
two categories s 

1) Ballistic Component: These electrons satisfy 

Equation (7) and have velocities l v ess than the velocity of escape 
such that 

2 

l/2mv - W - E^ < ecp Q 


and angular momentum require that 

"t = 

and 

l/2m(u ^+u ^)-ecp(r) = l/2m(.y r ^+v t ) -W 0 ~E£-ecp 0 (r) 


where R “r-s the radius of the body. 
Equations (4) and (5) yield 

2 2 ,, R?, 2 2 


u r - v +(1 - 7 2 ) v t “ {e(cp 0 -cn) + E f + W Q } 


( 8 ) 
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and hence describe ballistic orbits. 

2) Escape Component: These electrons satisfy Equation (7) 
and have velocities equal to or exceeding the velocity of escape 
such that 


l/2mv - W - E i em 
o • f o 


(9) 


and describe escape trajectories. 

These classifications are important in the evaluation of 
electron density and may be illustrated diagrammatically as 
in Figure 1. Curves I, II, and III describe equations 



v t i 2 - <|> CE f + W Q + ecp Q ] 


( 10 ) 


v r 2 + (1 - a 2 )v t 2 = (j|) f E f + W Q + efcp Q - cp}] , (11) 


and 

v r 2 = (H) [Ef + w o 3 


( 12 ) 


where a = R/r. These curves represent a circle, an ellipse 

and a straight line in the same order and distribute the thermionic 

electrons in various velocity doma'ins. 

The electrons with velocity domains external to the circle 

and the straight line - region A - are the escape electrons 

which do not return to, the body. The electrons with velocities 

* 

i 

i 

i 
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in the domain enclosed by the circle and the ellipse - region B - 

belong to the ballistic group with more than the necessary radial 

component of the velocity to reach the position r. These electrons 

are counted twice in calculating the electron density distribution 

1 

and make a dominant contribution to the local electron population. 

The electrons corresponding to region C - enclosed by the ellipse 
and the straight line - also belong to the ballistic group but 
do not possess enough radial velocity to reach position r. 

Therefore, these particles do not contribute to the local electron 
density. The straight line represents the least value of the 

v 

radial velocity that an electron must acquire before it can 
surmount the surface barrier. Therefore, the electrons corresponding 
to region D' in Figure 1 are not able to get out of the surface of 
the metallic body. 

III. ELECTRON DENSITY AS A FUNCTION OF POTENTIAL 

The contribution of the thermionic electrons with the 
initial (just inside the surface) velocities in the range 
(v, v + dv) to the electron population in a shell of radii r 

V 

and r + dr is determined by the product of the corresponding 
electron flux J(v)dv and the time dt = dr/u r spent by these 
electrons in traversing the thickness dr of the shell. This 
contribution d P t ^( r ) is given by 
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dp th (r) = a 2 f 


( 13 ) 


u r >0 u r 


which with the help of Equations (3) and (6) yields the 
expression for the thermionic component p^(r) °f the electron 
density at position r, 


p th (r ) = 4n (£) a2 I- 21 ! + ' 


( 14 ) 


where 


X. = 


j c 


v v dv dv 
r t r t 


[v r 2 +(l-a 2 )v t 2 -(2/m) fE f +W o +e(cp-co) }] 1/2 [e^ m(v r +v t ) " 2E f I I 2kT +l] 


Region B 


and 


(15) 


I = 
2 


v v dv dv 
r t r t 


J [v r 2 +(l-a 2 )v t 2 -(2/m) {E f +W o +e( 
Region A 


CD, 


,-«)}] 1 / 2 [e fm( V + V> - 2 E « 5 / 2k -hl ] 

( 16 ) 


and the limits of the integrals I, and I are set in accordance 

-he 2 

with Equations (10) - (12) and Figure 1, and the weight factors are 
inserted as explained in the preceding section. Including the 
contribution of the external plasma, the total electron density p (r) 
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at the position r becomes 


p(r) = p (r) + p (r) = 8n(mA) 3 a 2 (I, + I ) 
tn p 12 


(17) 


On introducing the following dimensionless parameters 


X 2 = mv r 2 /2kT ; Y 2 = mv t 2 /2kT ; e = E f AT 


(18) 


and 


a - [E + W + e(cp - cp) ]AT , 


Equation (17) reduces to 


p ( r ) = 8 n (2mkT) 3/2 g 2 f 

h 3 J 


XYdXdY 


[X 2 + (l-a 2 ) Y 2 -a] 1/, ^[e x ^ +Y ^” € +l] (19) 


But, 


X 2 +(l-a 2 )Y 2 -aaO 


2 2 

X +Y -e^a-es W./kT » 1 


and, therefore, we can neglect the unity term in comparison 
with the exponential term in the denominator of Equation (19) 

V 

Furthermore, on setting 


2 2 2 2 

X + (1-a ) Y - e - (W q At) = Z 


(20) 
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and 


(1 - a 2 ) x / 2 (y/x) = tan R 


(21) 


in the last equation, and after some simplification, the 
expression for the total electron density p (r) reduces to 

p (r) = 2 (2nmkTA 2 ) 3//2 e"^ W o +e (cn o -c P) 


[ 1 - fl-( R /r) 2 j 1/ V R2e <»o-®>AT(r 2 -R 2 )] 


( 22 ) 


which is expressed as a function of the potential cp(r) “ cp- 
y 

IV. REDUCED POISSON EQUATION 


Substituting for the electron density p(r) from Equation (22) 
in the Poisson Equation (1) and introducing the dimensionless 
quantities 

\tf = ecp/kT 

and (23) 

p = w At 

. o 


V 

we obtain 


d 


2 drV dry 


= A(kT) 1 / 2 e“^o“*^i-|l-(n/r) 2 ]. 1 / 2 e" (,|F o"* )R /(r " r2) ] 


( 24 ) 
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whore 

A = -2(e 2 /£ c ) (2 t®A 2 ) 3/2 e~ 9 (25) 

and the boundary conditions of the problem are 

*(R) = (26) 

at r = R, and 

iji (m) = 0 (27) 

at r = ®. 

In some cases of interest to us, we will find that the 

equilibrium potential energy ecp Q is much greater than the thermal 

energy kT corresponding to the surface temperature T. This would 

then enable us to neglect, to a first approximation, the second 
y 

term in Equation (24) , and hence we have 

7 2 \Jr (r) = Be X 

with 

1/2 

B = A (kT) , and X = i|f “ 

o 

Equation (27) is identical to the V so-called isothermal equation 
which has been solved for various boundary conditions and applied 
extensively to the problems pertaining to stellar structure by 
Chandrasekhar (1939) . , 


(28) 


(29) 


i 
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V. DETERMINATION OF SURFACE POTENTIAL 

The equilibrium value of the surface potential is determined 
from the balance of the escape component of the thermionic 
emission current and the plasma accretion current. The plasma 
accretion current consists of the electron and ion components. 

In the absence of a surface potential cf> o , the ion accretion 
current is smaller than the electron accretion current by a 
factor of the order of (m e /m^) . Therefore, only the relative 

initial magnitudes of the thermionic escape current and the 
plasma electron accretion current need be considered, and the 
ion accretion current may be neglected. Then, the surface 
potential cp 0 is positive if the initial thermionic escape current 
is greater than the initial electron accretion current. It may, 
however, become necessary to include ion accretion in consideration 
of the magnitude of the surface potential cn Q # if the latter is 
negative. 

Let us first consider the case of a positive surface 
potential. Then, the thermionic escape current is given by 

V 

v r v t dv r dv t 

J esc = 8ttR XX -fm(v r 2 + v 2) _ 2E f }/2kT . " . 

Region A I 
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ion IX, the electrons in region A of 


As mentioned in Section 

, a _j, transverse velocity components 
Figure 1 must have radial and transverse 

such that: 

V 2 > (2/m) (E + w ) 
r 1 


and 


v 


2 + V 2 > (2/m) (E f + W + ecp o ) 
r t 


Theref ore , the expression 
may he rewritten as 


for the thermionic escape component 


7(2/ taj (E f + W o + ecf 0 ) 


v t dV t 


J esc = 8ttr2 (n,A) / _ -{m(v7 + v7>- 2E £ }/2kT~ 


v 027mTT2 f +""vg /rz/STTSf 


(31) 


v dv 
t t 


+ X v r dv r £ -{mfv" 2 + v t 2 ) - 2 E f }/ 2 hT^ 


7(2/m) ' (E f + W Q +,ecp Q ) 


v 
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Once, again we can neglect the unity terms in comparison 
with the exponential terms in the denominators of the integrands 
in Equation (31) . After carrying out the integration and making 
some simplifications. Equation (31) reduces to 

J = (4iTEkT) 2 (m/h^) (l+(i) e ^ . (32) 


For bodies moving slower than the mean thermal speed of 
the plasma electrons, the electron accretion is symmetrical 
about the boc^y, and is given by 


J e * 


2rr 

3 


V»e v e R ' 


2 (T/T, 




(33) 


when n and v are the number density and the mean thermal 
e e 

speed of the plasma electrons, and r| e is the sticking coefficient 

defined as the fraction of the incident electrons transferring 

their charge to the body. In estimating the plasma accretion 

current J we further note that the ion accretion is further 
p 

reduced by a factor e and becomes negligibly small as 

compared to the electron accretion. Hence, 




2 (T/T ) # 
n n v,R P ° 

'e e e e 


( 34 ) 


which, when combined with Equation (32) in the condition of 
equilibrium, 
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J - J , (35) 

esc p 


yields 


[1 + t/t ]* 

L p Y o 

— i" 


24rm e (kT) e 


2 -0 


= 6 x 10 


29 T 2 e ~ 9 

r|n e v e 


(36) 


where T is expressed in electron volts. 

if * on the other hand, the surface potential cp Q is negative, 

' I (T/T ) iji I 

the ion accretion current is enhanced hy a factor of e p o . 

With | \jf 0 1 larger than a few tenths of an electron volt, the 

enhancement factor e 1 1 |/ 0 1 may be large enough to counteract the 

1/2 

effect of the reduction factor (m /m.) so that the ion current 

e i 

fnay by no means be negligible. In these circumstances, we must 
include the term 


J. = 

i 

in calculation of J 


( 2 rr/3 ) 


- (T/T 
e 


In writing Equation 


pHo 

(37) 


(37) 


we have assumed 


P ' 

that the ion accretion is also symmetrical about the body. If, ■ 
however, the speed of the body exceeds the mean thermal speed 
of the plasma ions by an order of magnitude, the ion accretion 
current (Equation (37)) is reduced by a factor of 1/2. The 
corresponding electron accretion current is given by 
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J e = (2tt/3) tl e n e v e R 2 e (T/ V . (38) 

i 

Therefore , the expression for the plasma accretion current , 


reduces to 



(2rr/3)R 2 n v [n e (T//a P^ ’" F o- n (m /m ) 
e e o i e j. 


l/2 e -(T/T p ) * 0] 
(39) 


In the calculation of the thermionic escape current we 
may first remark that the negative surface potential in our 
problem is only a fraction of a volt. It may also be noted that 
a negative surface potential, however small, enables all the 
emitted electrons to escape. Hence, the thermionic escape 
current is approximately given by 


esc 


2 3 - o 

= (4rrRh;T) (m/h ) e 


(40) 


Finally, in the condition of equilibrium (Equation (35) ) 
Equations (39) and (40) yield 


n a ‘e (T/ V^ -ri i (m e /m i ) 1/2 e (T/ V *0 = 24u (m/h 3 n e v e ) (kT) 2 e 6 


09 2 v -A 

- 6 x 10 (T /n v )e 
e e 


( 41 ) 
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VI. DISCUSSION 

k ' - 

In the preceding sections, we have formulated and analyzed 
the problem of the screening of the electric potential on a hot 
spherical object surrounded by an external plasma. It is assumed 
that 1) the spherical body acquires the electric potential in the 
processes of the thermionic emission of electrons ' from the surface 
of the object and the accretion of the charged particles from the 
surrounding plasma, and 2) the surface potential and the 
distribution of the potential and the electron density in the 
screening cloud are spherically symmetrical about the object. The 
basic requirement to satisfy these assumptions are that i) the 
surface of the spherical object is at a uniform temperature and 
ii) the object is either at rest or it moves with a speed that is 
small compared to the mean thermal speed of the plasma electrons. 
These requirements set restrictions on the exact application of 
the results of the present analysis to actual objects in space. 

The present analysis, nevertheless, provides, even in such cases 
where the above-mentioned assumptions do not strictly hold, at 
least an order of magnitude estimate of this phenomenon in front 
of the hottest part of the object. 

The applications of our analysis may be found in objects 


- 21 - 


on taring a planetary atmosphere or those approaching sufficiently 
close to a hot star. A space vehicle entering the earth s 

o 

atmosphere encounters stagnation temperatures of the order of 1500 K. 

o 

Ail meteoric objects acquire surface temperatures above 1200 K. 

Ionization in front of the cometary heads and certain cometary 

tails which is not understood as well, may be attributed in part 

to the solar heating of the metallic content of these objects. 

In general, the surface temperatures of the above-mentioned classes 

of objects are not uniform. Due to the variety of the types of 

such objects and uncertain available data, we will not make any 

attempt to apply our analysis to any specific case of the above- 

rnentioneo space objects. Instead, we will illustrate our theory 

by considering a hypothetical spherical object heated to a uniform 

surface temperature and surrounded by a plasma of electron density 

n — lO'Vcrn 3 at the equilibrium plasma temperature T ~ 1000 K (0.09 ev> 
e ■ & 

Two values of the work function and five values of the surface 
temperature, viz . , 

W '=3.0 and 3.8 (electron volts) 
o 

v 

and 

T = 0.04, 0.06, 0.09, 0.13 and 0.15 (electron volts) 
are considered to illustrate t he influence.of these parameters 
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on the nature of the electron cloud around a hot object. A 
common value of r ie = ^ = 0.1 is adopted for the sticking 
c ° 2 " ■ Llclen ^ :s ■ Since these surface-plasma parameters appear in 
a logarithmic term, any departure from this value for the 
SciC,\ing coerficients is not likely to seriously affect our 
result. 

T.v... equilibrium value - of the surface potential cr> Q is 

determined by the surface temperature T, the electron density 

n e' and temperature T D of the surrounding plasma. At low values 

0-L j. , thermionic emission of electrons is small, and hence, the 

balance of the electron and ion accretion currents from the 

surrounding plasma establishes a negative potential cp on the 

o 

objec v. s surfeice. Tne numerical value of cp Q is always a fraction 
(s 0.05) of a volt because even this small value of cp^ is large 
enough to increase substantially the ion accretion current and 

reauce the electron accretion current to off-set the relative 

_ _ 1 /•> 
errect of the factor (m /m ) . 

e i' 

Table 1 

SURFACE POTENTIAL OF A '•'METALLIC BODY 

W o ~ O.Sev T p = 0. 09ev (1044°K) n e = 10 3 /c.c rj - 0.1 

revolt) 

-0.1691 
-0.1683 
+0.0963 
+0.4395:' 

+0.8430 
+1.1340 


T (ev) 
0.04 
0.06 
0.09 
o.ir 
0.13 
0.15 
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At high values of T, on the other hand, a positive surface 
potential of several volts is established by the balance of 
the thermionic-emission and the electron-accretion currents; 
the icn-accre tion current having been reduced to a negligible 
value by the joint action of the positive potential and the 
factor (m /m ) 1//2 . Table' 1 lists values of cp corresponding 

e x o 

to the several values of W and T. 

. o 

Equation (24) can be reduced to a dimensionless 


differential equation, 


Q , ^ 

~2 dX X dX = ^ 


2 ;kt) x/ V < V*> Li- (i-x' 2 > 2 j 1/2 e - U 0 -*)/ix 2 -i) 

(42) 


where X = f . The variation of potential with distance from the 
spherical hot body of 1— ern radius is calculated by solving 
Equation (42), and the results are illustrated in Figures 2 and 3. 
The two curves representing the inclusion and the exclusion of 
the second term inside the parenthesis of Equation (42) for the set 
of parameters T = 0*09 ev, W = 3 ev and = 5.43 are shown in 
Figure 2, while Figure 3 exhibits the profile of the potential 
distribution, with the inclusion qf the second term in Equation (42) 
in the numerical calculations, for the two sets of parameters: 


i) T = 0.13 ev; 


W =3.8 ev ; and = 6.5, 
o o 
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and 

ii) T = 0.11 ev; W q = 3.8 ev; and $ o = 4.0 


respectively. The variation of potential with distance has 
the following characteristics : 

1) The nature of the profile of the potential distribution 
curve is independent of the set of the parameters used. The 

f 

potential falls very rapidly with distance from the object, 
and reduces to 1/3 of its surface value at a distance of 
approximately 2.3 and 1.7 cm in Figures 2 and 3 respectively. 

At a distance of about 8-10 cm the potential acquires an almost 
zero value and the surface potential of the body is completely 


shielded by an electron cloud of this dimension. 

2) The inclusion or the disregard of the second term inside 
the parenthesis of Equation (42) does not seem to matter in the 
calculation of the potential distribution. It is apparently due 
to the very rapid decrease of potential with distance from the 
object which reduces this term to a second order of exponential 
in - A thereby -making it negligible in comparison to the 

V 

first term. 


The elect 
body is caicu.L 
values of the 


ron density in the electron cloud surrounding the 
ated from Equation (22) by substituting in it the 
potential distribution obtained from the solution 
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of Equation (42) . The results of this computation are given in 
Figures 4 and 5. Figure 4, like Figure 2, includes two curves, 
one of these corresponds to the inclusion of the second term 
inside the parentheses of Equation (22) while the other disregards 
this term. The set of parameters used in the computation of 
these curves have the value T = 0.09 ev; W q = 3 ev; \ji Q = 5.43. 

The den. i tv distribution curve in Figure 5, as in Figure 3, 
corresponds to the sets of parameters having values 


i) T = 0.09 ev; W = 3 ev; i | i = 5.43 

o o 

and 

ii) T = 0.13 ev; W q = 0.33 ev; ijj 0 = 6.5 

and the numerical calculations are based on the inclusion of 
the second term in Equation (22). These curves bring out the 
following features of the variation of the electron density 
with distance from the object: 

1) There is a considerable increase of electron density 

in the immediate vicinity of the body. 

2) The electron density decreases very rapidly with 
distance from the body. 

3 ) Unlike in the estimates of the potential distribution, 
the inclusion or exclusion of the second term within the 
parentheses of Equation (22) in the computations of the electron 
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density appears to make a 'substantial difference in these 
estimates (see Figure 4) . The neglect of this term yields a 
value for the electron density at great distances which is 
higher than the ambient value. Therefore, it is necessary to 
consider this term in order to arrive at the correct estimates 
of the electron density. 

v;e wish to thank Dr.’ J. Herring for helpful discussions 
and Dr. H. Zapolsky for his helpful comments on some 
mathematical aspects of the analysis- 
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CAPTIONS TO FIGURES 


Figure 1 - Curves I, II and III are the plots of v r and in 

accordance with the Equations (12) , (13) and (14) respectively, 

and define .the velocity domains of the thermionic elections. 

Electrons in domain A form - the edcapd^grou?/' whereas those in 

coma: ' B and C describe ballistic orbits with sufficient and 

insufficient energies respectively to reach position r. 

Figure 2 - Plots of (j (r) against r/R. (W^ = 3ev; T = 0.09ev; 

and $ =5.43.) Solid line = 

o 

7 2 ^ = A (kT) [1- { 1- (R/r ) ** } exp(iJ(--(f 0 )F /(r -R ) ]exp ( 


2 1/2 

Dotted line = 7 ijr = A(kT) exp ( y ) 


Fxgure J 


3 - Solid line = W = 3.8ev; T = 0.13ev; A = 6.50, 


Dottea 


line = W = 3 . 8ev ; T = O.llev; * = 4.45. 


Figure 4 - (W = 3ev; T - 0.09ev; and * = 5.43.) Plots of 

o - ° 

log^ p(r) against r/R. Solid line = 

p = a[1-[1- (R/r) 2 } 1 /, 2 exp(’(f-^)R 2 /(r “R ) ]exp U" y Q ) • 


Dotted 3_.ine - p - A exp ( y * 


i* 


5 - Solid line = W = 3.8ev; T = 0.13 3v; and - 6.50. 


Dotted line = W = 3ev; T = 0.09ev; and ^ - 5.43. 









